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Abstract 

We examine various aspects of the normal state phase of Cu02 planes 
in the high Tc superconductors. In particular, within the context of the 
three band Hubbard model, we study as a function of doping the compe- 
tition between a charge density wave phase induced by oxygen breathing 
modes, antiferromagnetic order and paramagnetism. To account for the 
strong electronic interactions, we use the finite U slave boson method of 
Kotliar and Ruckenstein. 

I. INTRODUCTION 

Since the discovery of high temperature superconductivity in the cuprates, there 
has been a tremendous effort to deduce the properties of strongly interacting fermion 
systems in two-dimensions (2D). Regarding the copper oxide planes, this has led to a 
strong focus on purely electronic versions of the one and three band Hubbard models 
and their variants. Whilst this has led to progress in understanding the insulating 
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behaviour of the high Tc materials at half filling (as charge transfer insulators in 
the three band case and as Mott-Hubbard insulators in the one band case), there 
is little consensus as to whether superconductivity exists in either of these models or 
whether they are Fermi liquids or something more exotic. Furthermore, there have been 
relatively few attempts to understand the role of the underlying lattice in these systems. 
The motivation to study this last point becomes more acute with the realization that 
there is strong experimental evidence to suggest that a coupling does exist between 
charge carriers and the lattice in these materialsB^. 

From a theoretical viewpoint there have been several studies of both the one and 
three band Peierls Hubbard models in 2D. Here, we present a brief account of some 
recent results in this area. Dobry et al.i suggest using exact diagonalisation of a twelve 
site CU4O8 cluster, that in the stoichiometric state (zero doping), for any reasonable 
coupling of the fermions to the in-plane oxygen breathing modes and in the adiabatic 
limit, the system is stable against a lattice deformation. They also show that doping 
with one hole can produce an extended charge density wave (CDW) below a critical 
coupling to the lattice, and a self trapped polaron which forms a more generalised 
Zhang-Rice singlet above the critical coupling. On the other hand, in their Hartree 
Fock analysis of a system of 6x6 CUO2 unit cells, Yonemitsu et al.iil show that in the 
zero doped state the CUO2 planes can, above a critical coupling strength, undergo some 
type of lattice deformation, ie bond order wave, CDW. Slave boson studies of the one 
band Peierls Hubbard mo del0'0 show that at half filling a stable paramagnetic phase 
with an on-site frozen breathing mode induced long range CDW occurs below a critical 
value of the Hubbard repulsion U. Mitra and Behera showEi that in general electron 
correlation suppresses the CDW state. 

The purpose of this paper is to examine systematically as a function of doping the 
stability of the three band Hubbard mo dellim against a lattice deformation created 
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by the in-plane oxygen breathing modes. More to the point, we wish to study whether 
this model is susceptible to at least a paramagnetic CDW phase either at zero or finite 
doping, and more importantly, whether the copper oxide planes can develop from a 
homogenous antiferromagnetic phase in the undoped system to a paramagnetic CDW 
phase at finite doping. 

Using the slave boson method of Kotliar and Ruckenstein we are able to 

take to into account the strong on-site and nearest neighbour fermion interactions. 
This finite Hubbard U technique has the advantage of being able to cope with the 
strong coupling regime in a non-perturbative way, and at the same time it is in rea- 
sonable agreement with Monte Carlo resultJil even at the mean field level. Thus, 
armed with this apparatus, we are able to examine as a function of doping and in- 
teraction parameters at the mean field level, the competition between an undistorted 
antiferromagnetic state, the normal paramagnetic state and a paramagnetic CDW state 
induced by lattice deformations. Lastly, because we use the three band model, we are 
able to accommodate a more complete tight binding band structure which includes 
oxygen-oxygen hopping in order to produce a realistic Fermi surfacei^. 



II. MODEL CALCULATIONS 

To examine the possibility of a lattice instability, we include an electron-phonon 
coupling which modifies the copper-oxygen hopping strength depending on whether an 
oxygen ion is displaced towards a copper ion or away from iti^. To simplify calculations 
we make the following approximations. Firstly, we assume a frozen phonon picture for 
the oxygen ions and thus work in the adiabatic limit. Secondly, we assume that the 
oxygen-oxygen hopping is not modified by the relative displacements of the oxygen ions. 
With these provisions, the three band Hubbard Hamiltonian with electron-phonon 
coupling (three band Peierls Hubbard model) is given by 
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where c\„ creates a hole at a copper site i with spin cr and c]-^ creates a hole at an oxygen 
site i with spin a. The on-site copper (oxygen) energy is given by (cj). The hopping 
between copper and oxygen sites is represented by tij and its sign depends on the 
symmetry of the copper dj.2_y2 and oxygen and Py orbitals. The sum over (ij) refers 
to nearest neighbour copper and oxygen sites only. Without lattice deformation and 
in the nearest neighbour hopping picture, tij would be given by a constant t. However, 
tij is assumed to be modified by the in-plane breathing modes as stated above in the 
following manneriH. If the surrounding oxygen sites are deformed towards a copper 
site, then tij = t + auj, where Uj is the absolute value of the oxygen lattice displacement 
from its mean position. Similarly, if the surrounding oxygens are displaced away from 
a copper site, then the hopping becomes tij = t — auj. The parameter tjji in (|l|) 
represents oxygen-oxygen hopping and its sign also depends on the symmetry of the 
oxygen orbitals. The sum over {jj') refers to hopping between pairs of second nearest 
neighbour oxygen j and j' sites. The U term represents the Hubbard interaction at 
copper sites whilst V is the nearest neighbour interaction between holes at adjacent 
copper and oxygen sites. We neglect an on-site interaction at the oxygen sites since the 
probability of double occupancy of these sites by holes for physically reasonable hopping 
and interaction parameters and not too large doping is quite small. The classical force 
constant of the oxygen ions is represented by the parameter K. Following referenced!^ 
we introduce a dimensionless electron-phonon coupling constant A = a'^/Kt. 

To accommodate the strong on-site and nearest neighbour fermion interactions, 
we introduce four auxiliary bosons fieldsli^ for the copper sites. In this scheme, each 
possible electronic configuration of the copper site is represented by one of these fields. 
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The fermionic creation operator is mapped c}^ i-^ Q'io-cJo- where q^^ is defined as 



The boson fields Cj, Sia, and di represent empty, single with spin a and double occupa- 
tion of copper sites respectively. The square root terms in ensures that the mapping 
becomes trivial at the mean field level in the non-interacting limit U ^ 0,V 0. The 
KR slave boson method replaces the Hubbard interaction with a bilinear term involv- 
ing the boson double occupancy fields whilst renormalising the hopping energy of the 
fermion quasiparticles by the factor gj^. Unphysical states arising from the enlargened 
Hilbert space are eliminated via the following constraints; 

cLCia = Sl^ia + dldi (3) 

and 

eh + E 4.s,^ + did, -1 = 0, (4) 
(J 

which represent charge conservation and the completeness of the bosonic operators 
respectively. With this mapping in place the Hamiltonian becomes; 

{ij)a {jj')a i<7 

+ EscW + Ef^<H+ E n^U. + <H)c]vSv + E 

jo- i (jj>oo-' j 

- E 7.0(^1^.. + did,) + ^ 7:(ele, + ^ sls,^ + did, - 1) , (5) 

icr jo 

where and 7- are Lagrange multipliers enforcing the constraints on the boson fields. 

The partition function for the system can be expressed as a functional integral 
over complex (bosons) and Grassmann fields (fermions). Integrating exactly over the 
Grassmann fields leaves the partition function a 

Z = J D[e]D[s„]D[d]D[^]DW]e-''rr , (6) 



where the bosonic action is 

S':/;- = f^dr{ Tr In [t,^q,, + (7,. + + 9^)5,, 



rincr' 

+ E ^1(7: + 5. - 7.Js.. + E <Uf/ + 7: + 5. - 7.Jrf. 

+ E + - E 7: + E \ku] } (7) 

and the subscript nn imphes nearest neighbour summations only. 

To evaluate the partition function we employ the saddle point approximationlll 
where all the boson fields are c numbers. Accordingly, at T = the energy per cell is 
given by 

^ = ^ E ?^.^ii(cLs<. + ^-c-) + ^ E + h.c.) 

+ ^ E(e* + + ^ E S'^i- 

icr ja 

- ^ E l^Asl + + ^ E + E 4 + - 1) > (8) 

where N is the number of unit cells. 

Equation (^) must be minimized with respect to the boson fields. The CDW phase 
can be examined by dividing the copper lattice into two sub-lattices A and B, which 
results in a doubling of the unit cell and hence a halving of the Brillouin zone. The 
above energy per cell will be modified by inclusion of summations over the sub-lattices 
and by dividing the boson fields into e^,s^ and where v = A^B. Here we have 
assumed a paramagnetic CDW phase, thus s^a = s^a- 

To examine the purely antiferromagnetic state, the copper lattice is again divided 
and the unit cell doubled. The displacement of the oxygen ions is set to zero. Due to 
the symmetry of such a phase, = Sg^, jAa = 7sa- and = dg- The homogeneous 
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paramagnetic phase can be studied by a straight forward minimization of the above 
energy per celL 



A. Charge Density Wave Phase. 

In this subsection we derive the self consistent minimization equations for the CDW 
phase. Dividing the copper lattice into an A and B sub-lattices and minimizing the 
energy per CU2O4 cell with respect to Uj, we get 

1 d 
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(Aj)a (Bja) 



+ 4:Kuj^0 , (9) 



where M is the number of CU2O4 cells. Moving to Bloch functions the above equation 
becomes; 
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where 



^A^—y2( 
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= ^ E ( [e-^'^^^^^'^^^cUc,^,,) + h.c] - fe^'=^(^+"^)(cU.c,,,,) + h.c. 
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A. = ]^ E (- + h.c] + [^^.e-'=^(^'^^)(ci,,c,^,J + h.c. 

+ [-tkyC^'^^'^-^^Hcl^^c^,^,) + h.c] - [tkyC'^'^^'^-^^Hcl^^c^,^,) + h.c] ) (13) 



and 
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E ( [-^^.e-*^^(^+"^)(cUc,^k.) + h.c] - [iA:.e^'=^(^+"^)(cLk.c,,k.) + h.c 
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,fc^e^^.(|+«.)(4^^e^^^j+/,.e.]) . (14) 



The hopping parameters are defined by = t + auj and tg = t — auj. That is to say 
that the oxygen ions are breathed in towards the copper A sites thus strengthening 
the hopping between them, and breathed out from the copper B sites weakening the 
hopping between B sites and oxygen sites. The lattice displacement variable is in units 
of the Cu-Cu lattice spacing. The oxygen sites are labeled around a copper A site in 
the following way. Looking down upon the CUO2 plane, Xi is the oxygen site to the 
right of the A site, X2 is to the left, yi is above and 1/2 is below. 

Minimizing equation (|D with respect to the boson fields on the A and B sub-lattices, 
we obtain the following equations; 



tA^^ + Ud, + Wn^d^ - 27,d, + dX = , (15) 



+ SVnjS, - 27,s, + = (16) 



and 



t.^.^ + I'.e^ = , (17) 



where u = A,B. Combining these equations with the constraints we arrive at the 
minimization equations for the A and B sub-lattices; 



^(1-^) 



1 + 2dl - 2dl - 



+ 



:i + dl-n^)i^-dl) ^dli'f-dl)_ 
The total energy per cell is thus minimized by self consistently diagonalizing the fermion 
part of the effective Hamiltonian in order to calculate n^,,f2j,,A^, and uj and solving 
equations ([ToD and (|T^). 

In presenting results for this phase, we begin by studying the interplay between 
Coulomb interaction and electron-phonon coupling in the undoped state. We have 
used the parameter set0 of t = 1.3, tjjf = 0.65, V = 1.2, and ej — = 3.6 where all 
energies are in eV. The force constant K = 32t eV/A i. This parameter set is used 



for all calculations. Of course the overall stability of the paramagnetic CDW phase 
can only be considered after comparing the total energy to that of the homogenous 
paramagnetic and antiferromagnetic phases, which we will consider in later sections. 

Fig. 1 shows the charge density difference between copper A and B sites as a 
function of U for different electron phonon couphng strengths. In this CDW phase, 
the electron phonon coupling has increased the copper A site charge density due to the 
increased hopping amplitude t^- Conversely, the charge density at the B sites has been 
decreased significantly due to te- This process is hindered by the Hubbard repulsion at 
the copper sites, which naturally works against this increased charge density at the A 
sites. The results of Fig. 1 seem to indicate that at a certain value of U, for a given A, 
the system undergoes a first order transition to the homogenous paramagnetic state. 
At this point the Coulomb repulsion has become too strong to sustain the CDW. As U 
is increased form zero the A site charge density decreases, however the B site charge 
density does not increase at the same rate. There is an overall shift in charge to the 
oxygen sites as U is increased. This increase is shown in Fig. 2. The charge continues 
to build up on the oxygen sites until just before Ucrit, where there is a small reversal in 
charge flow. It can also be seen in Fig. 1 that the stability of the system as a function 
of U is quite sensitive to the electron-phonon couphng strength. Fig. 1 shows that for 
a reasonably strong coupling (A > 0.8) the system is quite robust (at least there is a 
local minimum in the energy) for any physically reasonable value of U for the Cu02 
planes. 

In Fig. 3 we show the difference in the charge densities on copper A and B sites 
A {= — rig) as a function of doping [5 = n — 1) for the parameter sets A = 1.00, 
U = 10 eV and A = 0.95, U — & eV. There is an almost hnear decrease in the 
CDW order parameter as the system is doped until one reaches 5 crit ~ 0.49 for both 
parameter sets. Prior to (5^^^^, the B site charge is almost localised due to an extremely 
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small effective hopping g^t^ ~ 0.1. The double occupancy on the B sites is such that 
(i| ~ 10~^, and remains at this magnitude as more charge builds up on the B sites 
as the system is doped. Because of this, the holes on the B sites become more and 
more localised as extra holes are added to the systems since — > as — 1 and 
dg As 6 ^ ^crit^ ^ 1) thus when 5 exceeds 5^^^^ and therefore > 1 
the system is required to obey the physical constraint d\ > — 1, causing a finite 
Hubbard energy on the B sites which makes the system in this phase unstable. 



B. Antiferromagnetic Phase. 

There have been previous slave boson calculations carried out on the antiferromag- 
netic and paramagnetic phases of the three band Hubbard model by Zhang et alii. 
However, they have used a parameter set containing a small oxygen-oxygen hopping 
parameter {tjji = 0.2 eV) and charge transfer energy (e^ — = 1.5 eV). Also, their 
value of U (6 eV), is just at the lower bound of what is a reasonable value for the 
Hubbard interaction on copper sites in the copper oxide planesH. This being the case, 
it is useful to redo and expand on some of these types of calculations with improved 
estimates of the hopping and interaction parameters, especially when comparing the 
energies of these phases to a phonon driven one. 

The saddle point equations in the antiferromagnetic phase are derived by again 
dividing the copper lattice into A and B sub-lattices. As stated before, the symmetry 
of this phase implies = 6^,5^1 = Sgi = Si;,Sai = = S| and dA = dg = d. 
Following a procedure similar to that in the CDW section, the following minimization 
equation is derived; 



U 



1 + 2(i^ — rii — rii^ 2d'^ — rii^ 



y(r^^n~Jn~^ |_y^(l + rf2 - ni){ni^ - d^) ^Jd'^{nn - d"^) 
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+ ^ 



1 + 2(P — Hi — rii^ 2(f — rii^ 



(19) 



where 



{Ai) {Bj) 



- i E((c\,c,.) + h.c.) = ±- UicUc,,) + h-c.) . 

(^i> (-Bi) 

and = = rig.n^-^ = = n^^ and = n^^ = rig^. Equation (|19D is solved self 
consistently with the diagonalized one body fermion Hamiltonian in order to minimize 
the energy per magnetic cell. 

Fig. 4 shows the quasiparticle band structure in the antiferromagnetic phase at 
zero doping using the reference parameter set and U = 10 eV. The coordinates in the 
figure represent the non-magnetic Brillioun zone. One can clearly see the lower and 
upper Hubbard bands and the dispersion of the oxygen bands in between. The lower 
Hubbard band is fully occupied at this hole concentration and consists of mainly (73%) 
of copper states. Doped holes go mainly to the oxygen sites and occupy states around 
the (7r/2,7r/2) region ( the magnetic Brillouin zone). It is interesting to compare this 
figure with the Hartree Fock band structure obtained by Yonemitsu et al.§. Although 
their parameter set is slightly different (in units of eV; t = l,tjj> = 0.5, U = 8, = 1 
and U oxygen = 3), the four lowest bands are similar with the same curvature and band 
minima. This is despite the fact that they include a Hubbard repulsion on the oxygen 
sites. For the same parameter set Fig. 5 shows the quasiparticle band structure for a 
doping of 5 = 0.25. It can be seen from this figure that new states have formed within 
the original insulating gap of the undoped system reducing the gap significantly. 

In Fig. 6 we show the dependence of the staggered magnetic moment m (= n^^ — 
n^J on U in the undoped state. As U decreases to ~ 1.75 eV the lower and upper 
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Hubbard bands begin to overlap. Beyond this point our results become inaccurate as 
the doping can not be pinned at zero. However, as U increases from this point there is 
a steady increase in the staggered moment until it saturates for large U at m ^ 0.75. 
Similar qualitative behaviour is found in the one band Hubbard model using Monte 
Carlo technique^. As shown by Zhang et a/.H the effect of tjj/ is to decrease the 
magnetic moment because it becomes energetically favourable for some of the holes to 
move on the oxygen network rather than occupy mainly copper sites. The behaviour 
of the staggered moment upon doping is shown in Fig. 7 for the reference parameter 
set. The dotted line represents U = 10 eV and the solid line is f/ = 6 eV. The local 
Cu moment vanishes almost discontinuously at Scrit indicating the likelihood of a first 
order phase transition. The precise point at which the transition occurs is unclear due 
to the extremely slow convergence near this region. It is worth noting that apparently 
U as little effect on the value of 6 at which the antiferromagnetic order vanishes. The 
variation of other parameters such as hopping and charge transfer energies have a much 
more significant effecH. The relatively small effect of varying U may be due to the fact 
that at least in the undoped state, the antiferromagnetic order parameter has already 
begun to saturate at f/ = 6 eV. 



C. Energy per Cu02 cell. 

In Figs. 8 and 9 we show the calculated energy per Cu02 cell as a function of doping 
for the CDW, antiferromagnetic and paramagnetic phases. The various energies in Fig. 
8 were calculated with an electron-phonon coupling of A = 1.00 and U = 10 eV. The 
figure clearly shows that for this parameter set the system is stable against a breathing 
mode induced paramagnetic CDW ground state and that as a function of doping the 
system moves from an antiferromagnetic state to a homogenous paramagnetic one. At 
this value of A the effective hopping parameter = 0.07. A stronger coupling drives 
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tg ^ causing the copper B sites to be completely localised and thus it becomes 
debatable whether it is useful to do calculations in such a region. However, if one 
believes that one can and assumes that the energy dependence on doping for the CDW 
phase has the same gradient as shown in Fig. 8 then it may be possible to find a 
stable CDW at Scrit > 0.4. The same type of calculations are shown in Fig. 9 but in 
this case U = 6 eV and A = .95. It is interesting to note that here the system starts 
in the antiferromagnetic phase at zero doping then moves to a paramagnetic phase 
at 5 ~ 0.33 and then to a paramagnetic CDW phase at 5 ~ 0.4. These calculations 
indicate that at least in the highly doped systems, the three band Hubbard model is 
susceptible at this mean field level to some sort of lattice distortion (we have shown 
here the paramagnetic CDW) given the appropriate interaction parameters. However 
we have shown that for a larger value of U and reasonably strong electron-phonon 
coupling, this model suppresses the CDW phase. 



III. CONCLUSIONS. 

Using the finite U slave boson method of krU to account for the strong fermion- 
fermion interactions, we have examined for different interaction parameters the possi- 
bility of an instability of the three band Hubbard model towards a paramagnetic CDW 
state induced by the oxygen breathing modes of the copper oxide planes. We have 
found that this model can at the mean field level, exhibit an instability towards such a 
phase at large doping given a low value of the Hubbard repulsion on copper sites f/(~ 6 
eV). For larger U the CDW is suppressed. Also, we have examined antiferromagnetic 
correlations at both zero and finite doping. We have shown that like the one band 
Hubbard model, the three band model antiferromagnetic order parameter saturates 
as a function of U and the system undergoes a phase transition to the homogenous 
paramagnetic phase at finite doping. A suggestion for further calculation would be to 
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examine antiferromagnetic correlations in the CDW phase, particularly its importance 
at low doping. 
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Figure Captions. 



Fig 1: Dependence of the CDW order parameter on the Hubbard repulsion U for 
different electron-phonon couphng strengths and zero doping. 

Fig 2: Total oxygen site occupancy for a CU2O4 cell as a function of U for different 
electron-phonon coupling strengths at zero doping. 

Fig 3: Behaviour of the CDW order parameter upon doping for different interaction 
parameter sets. 

Fig 4: Antiferromagnetic band structure for zero doping and t/ = 10 eV. The lowest 
band is fully occupied whilst the remaining bands are empty. 

Fig 5: Antiferromagnetic bandstructure at a doping of 5 = 0.25 and [/ = 10 eV. The 
original gap has narrowed and the second lowest band is partially filled. 
Fig 6: Behaviour of the staggered moment on Cu sites in the antiferromagnetic phase 
as a function of U at zero doping. 

Fig 7: Doping dependence of the staggered moment for different values of U . 

Fig 8: Energy per CUO2 cell as a function of doping for the various phases. Here 

[/ = 10 eV and A = 1.00. 

Fig 9: Same as Fig. 8 but here the interaction parameters are t/ = 6 eV and A = 0.95. 
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